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The Examiner is thanked for the indication that Claims 27 and 28 ate directed to 
allowable subject matter. The independent claims 23, 30. and.32 have been amended to include 
the subject matter of Qaim 27. and claims 27 and 28 have been cancelled. The remaining 
pending claims are believed to be allowable for at least the reasons that Claims 23, 30, and 32 
are allowable. 

The Office Action indicated one of the references listed on the PTO-SB-08 
(GRAMANN; ABF Algorithms Implemented at ARL:UT; ARL-TR-92.7; U. of Texas Appl. 
Res. Lab., May 1992) was not listed on the parent patent application (6.724.9 1 6), and therefore 
requested that Applicants provide a copy of the reference. A copy is enclosed as an attachment 
to this paper. However, please note that this document appears to be listed as the last reference 
on the first column of U.S. PatentNo. 6.724.916. Applicants request that the Examiner make 
this reference of record in the appHcation by initialing and returning the earlier submitted PTO- 
SB-08 listing this reference. No fee is believed to be due under 37 CFR 1.97 for this submission, 
however, should such a fee be required, kindly charge the fee or credit overpayment to Deposit ' 
Account 50-0281. 

In view of the above amendment, appbcant beKeves the pending application is in 
condition for allowance. An early indication of the allowability of the application in the form of 
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1.1 BACKGROUND 



Adaptive beamfonning (ABF) has been in use on several of the Heet sonar anays 
for many yeais. It has also proven to be a usefiil tool to the Sound Surveillance 
Underwater System (SOSUS) community, enhancing sonar performance through the 
reduction of overall beam noise, interference rejection, and consequenUy enhanced 
detection range. Tlius, ABFhas become an important asset in processing line aiiay data 
for passive detection problems. 



It is anticipated that ABF will play a significant role in future U.S. Navy acoustic 
systems. ConsequenUy. ARL:UT has implemented a variety of ABF algorithms. Several 
ABF algGrithms have been implemented on the ARL:UT AiUant minj-supeicomputer so that 
comparisons between these algorithms can be made on the same data, and the advantages 
(or disadvantages) of the different implementations explored. This document describes two 
of the ABF algorithms implemented, and some of the methods employed in these 
implementations. A previous documentl describes much of the processing performed after 
the ABF weights have been calculated. The reader is rcfened to this reference for those 
details. 
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2. ABF ALGORITHMS 
2.1 GENERAL DESCRIPTION 
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t^l^f^ ' ^^^"''^ ^ implementations. In a 

lo^T"' "f"''' ''''' "'"^ ^ ---ive 

c^M. ; 'r^^^'I'^'^^y **«P«°<1« only on the latest CSM. but on previous 

CSlv^ and the solutions associated with them. The cuirent Implementations allow a moie 
exact comparison between the algorithm solutions, and provide an estimate of the 
performance ofthe converged adaptive (i.e., realtime) system. Although most realtime 
systems use an adaptive implementation, the perf-pmmnce is expected to be close to that of 
the optimum solutions. ""ioior 

2.2 SINGLE POINT CONSTRAINT ABF 

TTie single point constraint adaptive beamformer (SP-ABF) is described and derived 
below. This beamformer is not expUciUy implemented on ARLrUT's AlHant computer 
Rather, its solution can be obtained ftom the beamformer implemented with a white noise 
gam constraint Hence, the Affiant implementation description is included in Section 2 3.2 
However, the SP-ABF is a basic minimum variance beamfomu^r, and its derivation shows 
Uie basic mathematical methods involved in the more compUcated cases. TTius its 
oenvanon is included here for completeness. 

2.2.1 Description and Derivation 

The SP-ABF algorithm is derived by representing the output of the array (i e 
beamformed output power) as 

PaBF^ =H'"/fH' , (2.1) 

^^TlslTclTL'^'''' (orcomplexconjuga.)^spose. 
subjec to a smgle constraint In this case, the look direction response of the array is 
const^d to unity gain. Restating the pxx,blemin mathematical term^ ' 

mm w"Rw subject to w"d = i , 2) 
where af is the look direction steering vector. 
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To solve this problem, the method of Lagrange multiplicxs is used. Forming the 
functional dm adjoins the constiaint to the function to be minimized gives 

F = w"Rw-X(w"d-l) , (2.3) 

Taking the derivative of the functional with respect to the weight vector and setting it to 
zero, one obtains 

~=/Jw-^»0 . (2.4) 

Rewriting gives 

w^XR~'d . (2.5) 
Since w^^d =1, one can solve for X, and write die fanal solution for die weight vector as 

R-^d 

This is die optimum beamf onner solution subject to die unity gain look direction constiaint 
This beamfonner is often referred to as a minimum variance distortionless response 
(MVDR) beamfonner. Unfortuhately, Uiis beamfonner sohition is very sensitive to sources 
of mismatch.2 Thus a more robust beamf onning solution is generally needed.3 The more 
robust beamfonner incorporates a white noise gain constraint, and is described in die next 
section. 

2^.2 Implementation on the AUiant 

The MVDR (widiout WNQ beamfonner has not been expUdtiy implemented on the 
Alliant Radier, if a user would like diis solution, die white noise gain constraint can be set 
to a very low value (e.g.. -100 dB). essentially disabling it. Additional details are provided 
below. 



5 



PACE • RCVD AT 1«/2»200$ 1 :41 :38 PM [Eastern Daylight Time] • SVR:USPTO-EFXRF-eO0 ' DNIS:2738300 * CSID:202 404 7380 ' DURATION (mm.ss):09-28 



10/28/05 FRI 13:45 FAX 20*2 404 7380 



NRL CODE 1008.2 PATENTS 



@)014 



^'^ constrStt^ constraint abf with white noise gain 

23.1 Description and Derivation 

"ITic SP-ABF with a white noise gain constraint is derived in a similar manner to the 
optimum beamform» described in the previous section. The idea behind this solution is to 
make the beamformer more robust, or tolerant of errors. One means of doing this is to 
constrain the white noise gain of the system. The white noise gain, for the MVDR 
beamfoimer, is defined as 



G 



1 



(2.7) 



The white noise gain is constrained to be a value greater than or equal to the white noise 
gain constraint (WNC), 52. By definition, the white' noise gain must be less than the 
number of elements in the array. Thus, when constrained, the white noise gain will faU in 
Che range 

S^^G^^N^ . (2.8) 

Thus the problem can be stated as 

nun w'^Rw aibject to w"d = 1 and S ff' . (2.9) 

The white noise gain constraint is a quadratic inequality consfraint, and thus the solution 
cannot be derived algebraicaUy in closed foim. The solution is the same as that of the 
previous optimum beamformer, except the CSM. represented by ^. is replaced by a CSM 
that has had "white noise" injected into it Thus, the CSM is represented by (R+el ), where 
£ ( the second Lagrange multipUer) is increased to the extent needed to satisfy the WNC. 
The solution now takes the form 
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Adding a small amount to the diagonal elements of the CSM (La.. £/) can be thought of 
m several ways. It de-emphasizes the contributions of the off-diagonal elements, which 
decreases the relative correlation seen between those coiresponding array elements. This 
makes the CSM appear to have weaker directional sources. Consequently the power 
nmumization does not attempt to "null" these contributions as much as it nonnaUy would 
without the increased diagonal contribution. Also, by increasing the diagonal contribution 
the CSM IS becoming more like a CSM of only white noise (diagonal matrix) As e 
approaches infinity, the CSM becomes diagonally dominant, and the resulting weights are 
the same as those for the conventional beamformer (CBF). This is the upper limit of the 
white noise gain. 

2.3.2 Implementation on the Alliant 

The implementation of the SP-ABF on the Alliant involves additional derivations 
and equations. First, the CSM can be written in terms of its eigenvalues (a) and 
corresponding eigenvectors (U ). Thus 

R^ULU" oiUdiag{a)U" . (2.11) 
The CSM with white noise injected can be written similarly, 

• R+el^U{l.^sI)U" , (2.12) 

with the inverse being 

(^+e^)-' = I/(Z+e/)-'y« . (2.13) 

nie eigenvalues and eigenvectors are calculated using a singular value decomposition 
routme (SVD). Rewriting the weights solution in temis of the eigemralues and eigenveaore 
gives . " 



UCT+£)j 



(2.14) 
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Rewriting the expression for 0^ in tenns of the eigenvalues and eigenvectors 



gives 



(2.15) 



This equation can be broken into the part that needs to be calculated once, and the part that 
varies wilJi the iteration on e Rewriting Eq. (2,15) gives 



(2.16) 



This expression (2.16) is first evaluated with csO. If it satisfies the WNC, then no fi is 
added to the CSM. If not, this expression is evaluated iteratively to find the proper value of 
e that satisfies the WNC. A version of Brent's method, adapted from Numerical 
Recipes, 4,* is used to find the value of s that satisfies the white noise constraint Since 
8 (as well as 0^), can span several orders of magnitude, e is parameterized in terms of 
log(fi) within the searching subroutine. Once a value of s is found that satisfies the 
constraint, the weights are evaluated using Bq. (2.14). The weights for one look direction 
are calculated for all frequencies and are written to disk, prior to calculating the next look 
direction. 

The white noise gain constraint is set relative to die maximum white noise gain 
possible for the array in use. In die current implementations, for frequencies above the 
design frequency of the array, the maximum white noise gain value is 10 Iog(Neiem). 
where Nelem is the number of sensors. This is the maximum gain attainable against 
spherically isotropic noise (Le.. the directivity index DI). At frequencies less dian the 
design frequency of the array, the maximum white noise gain is set to 10 log(2LaA). 
where La is defined as 



♦ Function ZBRENT, pp. 253-254. 



^"PaoEIbSFrCVD AT 10/28/2005 1:41:38 PM [Eastern Daylight Tlme|". SVR:USPTO-ePXRP-</30 ' DNIS:2738300 ■ CSID:202 404 7380 * DURATION (min^s):09-28 



10/28/05 FRI 13:48 FAX 202 404 7380 



NRL CODE 1008.2 PATENTS 



81017 



(2.17) 



and X are the sensor locations, L is the actual distance between ihe end sensors (i.e., Xmax 
and xmin), and X. is the wavelength at diat frequency. These define the maximum value of 
white noise gain, and the WNC is set relative to these values. Figure 2.1 shows an example 
of the white noise gain limit and a constraint curve -3 dB relative to this maximum. The 
airay has eight elements. Based on this constraint scheme, the white noise gain will fall on 
or between the two curves. 
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Figure 2.1 shows the current white noise gain constraint implementation. This 
implementation is a compromise solution that is cuirently being used in the existing 
software. Several conections to this constraint (e.g., setting die maximum white noise gain 



9 



PAGE 17/29 * RCVD AT 10/28/2005 1 :41 :38 PM [Eastern Daylight Time] • SVR:USPTaEFXRF4/30 * DMS:2738300 * CSID:202 404 7380 * DURATION (tnm-SS):09-28 



10/^8/05 FRI 13:47 FAX 202 404 7380 NRL CODE 1008.2 PATENTS BlOlS 



limit to 0 dB for very large wavelengths, accounting for variations in steering direction, 
etc.) will be examined in fixture work, and reported in a later document 

2.4 MULTIPLE LINEAR CONSTRAINTS WITH WHITE NOISE GAIN 
CONSTRAINT 

2,4,1 Description and Derivation 



Multiple linear constraints can be added to the beamformcr solution to provide 
additional control of the beamformer response. These additional constraints can be used to 
force the beam response curve to have its first, and possibly higher derivatives equal to 
zero in the look direction. This type of linear constraint is known as a derivative constraint 
Thus one can set the first derivative to zero in the look direction, forcing the beam response 
curve (i.e., the beam pattern) to be relatively flat in tiie vicini^ of the look direction. 

Another type of linear constraint is the point constraint This constraint is used to 
set a specified response at locations other than the look direction. Point constraints are 
generally used as maiidobe maintenance constraints. One example is the "three-point 
constraint," where the look direction and one point on either side of the look direction is set 
to a specific level. Another constraint scheme mi^t be to match the CBF beam response at 
specific locations on either side of die look direction. This is referred to as the Match-CBF 
scheme. The user selects the points (in dB) at which they wish to have the ABF beam 
response match the CBF beam response. This defines the location (in sin B space) that 
these constraints are applied. They are calculated for the specific array and the locations of 
the constraints vary with frequency. The current implementation on the AlUant allows the 
user to select the response points to match the CBF pattern, and then iieratively finds the 
locations of those constraint points. A "three point" scheme could be to match the CBF 
response where the CBF response is "3 dB down." Another possible "three point" scheme 
could be to match the CBF response at the main response axis (MRA) of the adjacent 
beams in a multi-beam system. This scheme has not been implemented to date, but could 
easily be incorporated in the existing code. One can also specify response levels at locations 
relative to the look direction, which are independent of frequency. These are ^cUy a "user 
defined" point constraint 
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The multiple linear constraint APp {MLC-/:BF) is derived in the same fashion as 
the SP-ABF. In tiiis case, however, the look direction constraint iy^d =1) is replaced with 
a constraint niatrix and a response vector. Thus the single constraint becomes 

= (Z18) 

where Af is the constraint "location" matiix. w is the weights vector, and ^ is the response 
vector. 

Derivative Constraints 

M is constructed differently for derivative constraints than for point constraints. 
For derivative constraints, the first column of Af is the steering vector for the look 
direction. Suhsequait columns coirespond to the derivatives of the steering vector (i.e„ the 
second column is the 1st derivative, the third is the 2nd derivative, etc.), up to the number 
of derivatives being set to zero. Thus, the first column is 

where /is the frequency, c is the speed of sound, and 0 is the look direction. The second 
column contains the 1st derivative, and is calculated using the following expression: 

A/i.2 = exp^-y^^)M,.t . (220) 

Higher order derivative columns then become 

^k.«d^rh-^l='Mt,2Mt^^ . (2.21) 

The response vector associated with this consttaint matrix has 1 for its first element, 
corresponding to the look direction response, and zeroes for all of the derivative terms. 
Thus,; takes the form 
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8^ 



1 
0 
0 



(2.22) 



with as many zeroes as derivatives being set. 



Point Constraints 



The point constraint matrix is constructed with the columns forming the direaion 
vectors for the constraint locations- Thus, the first column is the look direction, as before, 
with subsequent columns forming the looJc direction vectors for the additional point 
constraints. Hence, 



(2.23) 



where 9n are the different constraint locations. The response vector, g, contains the look 
direction constraint (i.e., 1) in the first element, with the other constraint responses 
following. For instance, if three constraints were set (i.e., look direction plus two point 
constraints) and these additional point constraints were set to 3 dB down, g would be 



LOO 

0.707 

a707_ 



(2.24) 



Solution for Multiple Linear Constraints 



The derivation of the solution is similar to that for the SP-ABF. Stated in 
mathematical terms, the problem is 



min w^Rw subject to M"w^g and G^^S^ 



(2.25) 



Forming the functional with the multiple linear equality constraints gives 
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F = w«Rw-X{M»w-g) . (2.26) 

Taking the derivative with respect to the weight vector, and solving for A, gives the 
solution for the weight vector, 

w = R-'MlM"ir'Ml'g . (2.27) 

This equation shows the solution without the white noise injected. Injecting white noise 
replaces R with {R +d ). Thos, the sohition with the WNC included is 

w = iR + El)-'M[Ad''(R + Bir'M]~'g . (2.28) 

2.4.2 ImpIemeatatioD on the Alliant 

The MLC-ABF has been implemented on the AUiant for the three types of constraint 
schemes described above (i.e., derivative, match-CBF, and point constraints). Only one 
type of constraint can be used at a time with the cuirent code. As witii flie SP-ABF case, 
the CSM is decomposed into eigenvalues and eigenvectors using an SVD subroutine 
(Eq. (2.11)). Based on this decomposition and the solution for the weights, the expression 
for the inverse of the white noise gain (i.e., I/Gw) can be written explicitly as 

y^''yv^[{Udiag[a+Er'U")M{M»Udiag[a+ epU'^Mf's^" x 

{{Udiag[a+E]-'u'')M{M«Udiag[a+epU''M)''g^ . (2.29) 

This equation can be rewritten as 

>v*vv = {g''lM"Udiag{a+e)-'u''M]~'M''u"diagia+ e)-»t/| x 

{udiagi<r+e)-'u"Ad[M"Udiag(a+er'u"M]''gj . (2.3O) 

This equation is broken into components that can be computed once, and those changing 
with iteration. Two terms computed once and stored are 
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A = V^M and = {M^J]\fyU^M) 



Consequently, 



(2.32) 



is computed as 



1 



(2.33) 




Again, the iteration on white noise gain is perfonncd if the solution without white noise 
injection does not satisfy the white noise gain constraint Prior to iteration, the bounds of e, 
which is parameterized as log(£), are found using an initial bracketing routine adapted from 
a subroutine in Numerical Recipes A** After the initial bounds arc set, a root bracketing 
function using Brenfs method (i.e,, the same as SP-ABF) is used to find the proper £ 
value, and the weights are calculated and stored for that frequency. After completing all 
frequencies in the CSM, the weights for that look direction are written to disk. 

2.5 EXAMPLE BEAM PATTERNS AND SOLUTIONS 

Figures 2.2 through 2.4 show examples of beam patterns calculated from the ABF 
weights for the different algorithms described above. All were calculated using the same 
four CSMs and the same array configurations. The CSMs were cakulated from timeseries 
data obtained from a horizontal line array in a deep ocean environment The time period 
between each CSM is approximately 15 s. In all cases, the white noise gain constraint was 
turned off (WNC -100 dB). The range of sin 0 in each figure is -2 to 2, thus showing the 
beam pattern in visible and invisible space. TTiis allows inspection of the sidelobe stmcture. 
Also, mainlobe squinting at look directions near endfire can be seen more easily with this 
format. 



* Subroutine ZBRAC, pp. 245-246. 
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Figure 2.2 shovrs the beam patterns based on the SP-ABF, and the equivalent CBF 
beam pattern. The CBF beam pattern is offset -10 dB for clarity. aeaily» each of the beam 
patterns satisfies die look direction constraint of unity gain (0 dB response). However, the 
slope of the response curves in the look direction, as well as the shapes of the mainlobes 
vary with the different input CSMs. 

Figure 2.3 shows the beam pattern calculated from weights generated with the 
MLC-ABF program, with the first derivative set to zero in the look direction. The slope of 
each curve is zero at the look direction, and the response xs much flatter (and closer to 
0 dB) at locations near the look direction than in Figure 22. 

Finally, Figure 2.4 shows flie beam pattern calculated with weights generated with 
the MLC-ABF program using the Match-CBF option. In this case, a "three point*' 
constraint was set with the two mainlobe maintenance constraints set to match the CBF 
response at the "3 dB dovm" points. Now, all of the response curves have essentially the 
same mainlobe shape between sin 9 = -0.4 and 0.4. 
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3. SUMMARY AND CONCLUSIONS 

Two adaptive beamfonning (ABF) algorithms have been implemented at ARLrUT 
in order to assess the achievable performance of ABF, The single point (look direction) 
constrained (SP-ABF) and multiple linear constrained (MLC-ABF) adaptive beamformeis 
were derived and explanations of their implementation were given in the previous sections. 
A discussion of the white noise gain constraint and its implementation was also included. 
Additional wozk examining the performance of these algorithms is under way using data 
obtained with ARL:UT acoustic data acquisirion systems, and will be reported in a later 
document 
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